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Abstract 

About 2.9 fbr 1 data set had been collected with the BESIII detector at the BEPCII collider. Taking advantages of 
data near open-charm threshold, we present preliminary results on leptonic and semileptonic charm-meson decays, as 
well as result on the FCNC process D° — > yy from BESIII experiment. High precision charm data will enable us to 
validate Lattice QCD calculations at the few percent level. These can then be used to make precise measurements of 
CKM elements, V c d and V cs , which are useful to improve the accuracy of test of the CKM unitary. 
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1. Introduction 

The charm physics potential at open-charm thresh- 
old plays important role on quark flavor physics (T). 
The open-charm program at BESIII includes studies 
of leptonic, semileptonic, hadronic charm-meson de- 
cays, searches for neutral D mixing, CP violations, and 
rare and forbidden charm decays, which are sensitive to 
physics beyond standard Model (SM). The quark mix- 
ing parameters are fundamental constants of the SM. 
They determine the nine weak-current quark coupling 
elements of the CKM matrix [2 1. Studies of the semilep- 
tonic or pure-leptonic decays of D mesons are preferred 
way to determine the CKM elements, \V CS \ and \V c d\, 
since the strong interaction binding effects are param- 
eterized by form factors or decay constants that are cal- 
culable, for example, by lattice QCD (LQCD) and QCD 
sum rules. On the other hand, |V* CS ] and \V c j\ are tightly 
constrained when CKM matrix is assumed to be uni- 
tary. Therefore, measurements of charm semileptonic 
or pure-leptonic decay rates rigorously test theoretical 
prediction of the D meson semileptonic form factors or 
decay constants. High precision predictions of QCD 
will then remove road blocks for many weak and flavor 
physics measurements, such as in B decays for determi- 
nations of \Vub\ and neutral B mixing parameters. 

Many of the measurements related to charm decays 



are also accessible to the B-factories and Super-B fac- 
tories. What are the advantages to running at the open 
charm threshold at the BESIII experiment? The BE- 
SIII experiment will not be able to compete Super-B 
factories in statistics on charm physics, especially on 
the rare and forbidden decays of charm mesons. How- 
ever, data taken at charm threshold still have power- 
ful advantages over the data at T(45), which we list 
here JT] [3J: 1) Charm events produced at threshold 
are extremely clean; 2) the measurements of absolute 
branching fraction can be made by using double tag 
events, which was first applied by the MARKIII Col- 
laboration at SPEAR H. The produced ^(3770) in our 
sample decays into a pair of DD. Reconstructing one 
of the D mesons with know exclusive hadronic modes 
while looking for decays of the other D mesons would 
allow us to reconstruct final states with neutrinos; 3) sig- 
nal/Background is optimum at threshold; 5) Quantum 
coherence allow simple and complex [6| methods to 
measure D°D° mixing parameters, direct CP violation 
and strong phase differences. 

In this talk, I report preliminary results on charm me- 
son decays based on a sample collected at the BEPCII 
with the BESIII detector |7|. This sample was collected 
at s[s — 3.773 GeV with an integrated luminosity of 2.9 
fb _1 in which the background levels at the open charm 
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threshold is expected to be substantially lower than that 
at the T(45) peak at B factories. 

2. Purely leptonic D decay 

With a sample of 2.9 fbr 1 taken at open-charm thresh- 
old, BESIII experiment measures the decay branch- 
ing fraction for D + — > /i + v /J and extracts decay con- 
stant fo+ JH. Decay constant characterizes the strong- 
interaction physics at the quark-annihilation vertex. In a 
fully leptonic decay, they parameterize all of our essen- 
tial theoretical limitations. Decay constant for B case 
also appears in the evaluation of box diagrams, and limit 
theoretical precision in calculating the neutral B meson 
mixing. Thus, lack of knowledge of the B° and B s de- 
cays constants limits the usefulness of precise measure- 
ments of B° - B° and B s - B s oscillations. These mixing 
data are our best sources of information on the CKM 
matrix elements Vtd and V ts , which are difficult to mea- 
sure directly in top decay. The leptonic decay of charm 
meson presents an opportunity to check LQCD results 
for decay constants against precision measurements. 




1.82 1.84 1.86 1.88 1.82 1.84 1.86 1.88 1.82 1.84 1.86 1.88 



M B [GeV/c 2 ] 

Figure 1 : Distributions of the beam energy constraint masses for D~ 
hadronic tags used for the D + — » p + v analysis. Modes in panels (a)- 
(i) are D~ -> K^tCiC, K s iC \ K S K~, K+R~n~, TP/Cn'/P, t&tCiC, 
K s n~n a , K + n~n~n~n + ', and K s n~n~n + . 

A tag is simply a fully-reconstructed D hadronic de- 
cays. A sample of tagged events has greatly reduced 
background and constrained kinematics, both of which 
aid studies of how the other D in the event decays. One 
can infer neutrinos from energy and momentum con- 
servation, allowing full reconstruction of (semi)leptonic 
D decays. The typical tag rates per D (not per pair) 



are roughly 15% and 10% for D° and D + , respec- 
tively. For pure leptonic decay in this analysis, the 
singly tagged D mesons are reconstructed in nine non- 
leptonic decay modes of D — > K + n~n~, K s tC, K s K~, 

K + K~7C, K+TCTCvP, T^Tnr, K s n~n°, /T 7T- K if "7T + , 

and K s jt~j:~ji + . Mass peaks for the nine hadronic tag 
modes are shown in Fig. [T] A maximum likelihood fit 
to the mass spectrum yields the number of the singly 
tagged D events for each of the nine modes. The total 
number of tagged D events are 1565953 + 2327 D . 

The chosen signal variable for the jj. + v decay is the 
calculated square of the missing-mass of any undetected 
decay products, shown in Fig. [2] this should of course 
peak at M„ = for signal events. The power of D- 
tagging is evident in the clean, isolated signal peak. In 
Table [T] sources of background modes are summarized 
including Kin + , n + jr°, rv and other components. After 



Table 1: Sources of background events for D + — > pv. 



Source mode 


Number of events 


D + -» K L n + 


7.9 + 0.8 


D + -> 7T+7T 


3.8 ± 0.5 


D + -> T + V 


6.9 ± 0.7 


Other D decays 


17.9 ±0.1 


ISR to f and J /if/ 


0.2 ± 0.2 


i/<3770) -> Non - DD 


0.9 + 0.4 


e + e~ — > qq 


8.2+ 1.4 


e + e > t + t~ 


1.9 ±0.5 


Total 


47.7 ± 2.3 
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Figure 2: BESIII missing-mass-squared plot for D + — » p + v. The in- 
sert shows the signal region on a vertical log scale, where dots with er- 
ror bars are for the data, histograms are for the simulated backgrounds 
from D + -> K L n + (red), D + -> n + n° (green), D + -> r+y (blue) and 
other decays of D mesons (yellow) as well as from e + e~ — > non-DD 
decays (pink). 

subtracting the expected number of background events, 
about 377.3±20.6 signal events for D + — > fi + v decay are 
retained, where the error is statistical. The overall effi- 
ciency for observing the decay D + — > /u + v is determined 
to be 63.82% by analyzing full Monte Carlo (MC) sim- 
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ulation events of D + — > fi + v versus D tags. Therefore, 
we obtain the preliminary branching fraction to be 

&(D + -> n + v) = (3.74 + 0.21 ± 0.06) x 10~ 4 , (1) 

where the first error is statistical and the second sys- 
tematic. This measured branching fraction is consistent 
within error with world average of S(D + — > fi + v) - 
(3.82 ± 0.33) x 10~ 4 [9], but with more precision. The 
decay constant fa is then obtained using 1040 + 7 fs as 
the D + lifetime and 0.2256 as V cc i [9|. Our preliminary 
result is 

fa = (203.91 + 5.72 ± 1.91) MeV, (2) 

where the first errors are statistical and the second sys- 
tematic arising mainly from the uncertainties in the 
measured branching fraction (1.7%), the CKM matrix 
element V a t (0.3%), and the lifetime of the D + meson 
(0.7%) $9). The total systematic error is 1.0%. In Ta- 
ble^ the BESIII preliminary results are compared with 
CLEO-c ED and LQCD predictions Ol for the decay 
constant. The preliminary result from BESIII is the best 
number in the context of the SM, and remains consistent 
with LQCD. 



Table 2: Comparison of results for D + — » fiv and decay constant from 
BESIII 1 1 1, CLEO-c 1 10 ) and LQCD prediction (TT). 



Model 




S(D + -^p + v)x 1(T 4 


f D+ (MeV) 


BESIII 




3.74 ±0.21 ±0.06 


203.91 ± 5.72 ± 1.97 


CLEO-c 


ED 


3.82 ±0.32 ±0.09 


205.8 ± 8.5 ± 2.5 


Average 




3.76 ±0.18 


204.5 ± 5.0 


HPQCD 


mi 




213±4 



3. Semileptonic D decays: D° -» K e + v and n e + v 

One of the best ways to measure magnitudes of CKM 
elements is to use semileptonic decays since they are 
far simpler to understand than hadronic decays and the 
decay width is ~ \V C<] \ 2 . On the other hand, measure- 
ments using other techniques have obtained useful val- 
ues for V cs and V c d fl2l . and thus semileptonic D decay 
measurements are a good laboratory for testing theories 
of QCD. For a D meson decaying into a single hadron 
(h), the decay rate can be written exactly in terms of the 
four-momentum transfer defined as: 

q 2 = (f D - ]f h f = m 2 D + m 2 h - 2E h m D . (3) 

For decays to pseudoscalar mesons and virtually mass- 
less leptons, the decay width is given by: 



where pp is the three -momentum of pseudoscalar me- 
son in the D rest frame, and f+(q 2 ) is a form-factor, 
whose normalization must be calculated theoretically, 
although its shape can be measured. 

The BESIII experiment has taken about 2.9 fbr 1 data 
at open-charm threshold during the 2010 and 201 1 data 
runs. Using one-third of the data, a partially-blind anal- 
ysis has been done with the D° — > Kev and D° — > 
nev decays. Using the double tag technique, several 
hadronic D decays are fully reconstructed at first. The 
following four hadronic D decays are used: D° — > 
K-jt + , K-jt + 7t 1 \ K-7T + 7r°7T° and K-^n-^. Mass peaks 
for the four hadronic tag modes are shown in Fig. [3] 
After hadronic D° tags are found, we reconstruct signal 




Figure 3: Distributions of the beam energy constraint masses for D 
hadronic tags used for the D° — > Kev and nev analyses. Modes in 
panels are D° -> R-n + , R-n+n , K-n + n {) n° and R-n+n-n*. 

decay for the other D°. The signal events with a missing 
v are inferred using the variable U = E m j SS - \P m i SS \, sim- 
ilar to missing mass square, where "miss" here refers 
to the missing energy or momentum. Figure. [4] shows 
the U distributions and fit projections for the decays of 
D° K + e-v and D° -> Tr+e^v. 

Given the signal yields obtained from fitting U dis- 
tributions and signal efficiencies obtained from signal 
Monte Carlo, the absolute branching fractions are ob- 
tained. Preliminary results of branching fractions are 
listed in Table [3] and comparisons with results from 
PDG2012 [9| and CLEO-c results 031 are also made. 
In order to measure form factor, partial decay rates are 
measured in different q 2 bins. D° — > K + e~v candidates 
are divided into nine q 2 bins, while D° — > n + e~v can- 
didates are divided into seven q 2 bins. Signal yields in 
each q 2 bin are obtained by fitting U distributions in that 
q 2 range. Using an efficiency matrix versus q 2 , obtained 
from Monte-Carlo simulation, and combining with tag 
yields and tag efficiencies, the partial decay rates are ob- 
tained, as shown in Fig. [5] The values of q 2 -dependent 
form factors in each q 2 bin can be extracted from the 
measured partial decay rates as shown in Fig. [6] These 
data can be fitted with different parameterizations of the 
form factors, and the fit can distinguish between form 
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Figure 5: Partial decay rates measurement using individual tag modes (points) and all tag modes combined (histogram) for decay of D° — > K + e v 
(left) and 5° -> n + e~v (right). 



Table 3: Branching fraction measurement using 923 pb 1 of 1^(3770) 
data from BESIII experiment, and comparisons with results from 

CLEO-c 1 13 | and PDG20 1 2 1 9 1 . 

Experiment S(D U K + e~v)(%) S(D" -» n + e~v)(%) 
BESIII 3.542 ±0.030 ±0.067 0.288 ± 0.008 ± 0.005 

CLEO-c [13] 3.50 ± 0.03 ± 0.04 0.288 ± 0.008 ± 0.003 

PDG2()12|9| 3.55 ±0.04 0.289 ± 0.008 



factors parameterizations. In general, one may express 
the form factors in terms of a dispersion relation, an ap- 
proach that has been well established in the literature 
(see, for example, Ref. [ 14 1 and references therein): 

- r M ^ dt , ( 5) 

n J(m D +m P y- t-q l -le 

where m po i e is the mass of the lowest lying (qiq/) meson 
with the appropriate quantum numbers: for D — > Kev 
it is D* + and for D — > nev it is D* + , the parameter a 
gives the relative contribution from the vector pole at 
q 2 = 0, mo is the mass of the D meson, and mp is 
the mass of the final state pseudoscalar meson. The 
integral term can be expressed in terms of an infinite 
series [14|. Typically it takes only a few terms to de- 



scribe the data. Three different parameterizations of the 
form factor f+(q 2 ) are considered. The first parameteri- 
zation, known as the simple pole model, is dominated 
by a single pole lfl5ll : the second parameterization is 
known as the modified pole model lfl5ll : the third pa- 
rameterization is known as the series expansion |[T4l . 
Thus minimized x 1 fits are employed to extract the val- 
ues of f+(Q)\V c d(s)\ using each of the parameterizations. 
The preliminary results for /+(0)| V C£ /( S )| are shown in Ta- 
ble^ With = 0.2252 (| V„| = 0.97345) and BE- 
SIII new results (3 par. series) [16|, we extract the val- 
ues for f+^ n (0) and f^~* K (P), and results are compared 
with other experiments and theoretical calculations as 
shown in Fig.[7][17|. 

4. Rare charm decays at BESIII: preliminary re- 
sults on D° -» yy 

Searches for rare-decay processes have played an 
important role in the development of the SM. Short- 
distance flavor-changing neutral current (FCNC) pro- 
cesses in charm decays are much more highly sup- 
pressed by the GIM mechanism than the correspond- 
ing down-type quark decays because of the large top 
quark mass. Observation of FCNC decays D — > hl + l~ 
and D — » Z + T could therefore provide indication of new 
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Figure 6: Plot of f+(q 2 ) from data with Lattice QCD prediction. Black 
points are from data with statistical error only, red curve is theoretical 
prediction, and blue curves are theoretical variations within statistical 
error only. The top plot is for 5° — > K + e~v and the bottom plot is for 
D° -> 7r + fTv (bottom). 



Table 4: Results of f+(0)\V c j( s )\ from individual form factor fits; sta- 
tistical and systematic uncertainties on the least significant digits are 
shown in parentheses. Results from CLEO-c [131 are compared. 







BESIII 


CLEO-c 


3 par. Series 






D° -> K + e-v 


0.729(8)(7) 


0.726(8)(4) 


D° -> 7T + e-V 


0.144(5)(2) 


0.152(5)(1) 


2 par. Series 






D° -> K + e-v 


0.726(6)(7) 


0.717(6)(4) 


D° -> n + e-v 


0.140(4)(2) 


0.145(4)(1) 


Modified pole 






D° -> K + e-v 


0.725(6)(7) 


0.716(6)(4) 


D° -> n + e-v 


0.140(3)(2) 


0.145(4)(1) 


Simple pole 






D° -> K + e-v 


0.729(5)(7) 


0.720(5)(4) 


D° -> 7i + e-v 


0.142(3)(1) 


0.146(3)(1) 



Figure 7: Comparison of f+(0) from experiments and theory. Left plot 
is for /f -"(0), and right plot is for /?"* x (0). Note that BESIII result 
from D decays only, while CLEO-c use both D° and D + decays 1171 . 



physics or of unexpectedly large rates for long-distance 
SM processes like D — > hV, V — > l + l~, with real or 
virtual vector meson V. 

With 2.9 AT 1 data at (A(3770) peak, BESIII search 
for D° — > yy decay which must be produced by FCNC. 
From the short distance contributions, the decay rate for 
D° ^ yy is predicted to be 3 x 1CT 11 US] [13 HDJ. 
However, the long distance contributions significantly 
enhance the decay rate which is estimated to be (1 — 
3) x 1(T 8 US HO). This decay could be enhanced by 
new physics (NP) effects which lead to contributions at 
loop level [21,22]. For instance, in the framework of the 
Minimal Supersymmetric Standard Model (MSSM), the 
calculation shows that the decay rate for c — > uy tran- 
sition could be 6 x 10~ 6 , which is one to two orders of 
magnitudes enhanced relative to the SM rate, by consid- 
ering gluino exchange ED . 

Experimental searches for D° — > yy were performed 
by the CLEO EJ] and BABAR (24) experiments based 
on data samples collected at the T(4S) peak. They 
found no significant signals. The latter experiment 
yields the most stringent experimental upper limit to 
date on the S(D° -> yy), 2.2 x 1(T 6 at 90% confidence 
level (C.L.). 

Since one of the major backgrounds in the analysis 
of D° — > yy comes from D° — > 7r°7r°, we firstly mea- 
sure the decay rate of D° — > 7T°n° using the same data 
sample. Then we use the measured D° — > tt 7t° decay 
rate to do the normalization and background estimation 
for D° — > yy measurement. Figure [8] shows the beam- 
constraint mass distribution of the observed 7t°7r sig- 
nal events and comparison with MC simulations of the 
expected backgrounds. The fit to the beam-constraint 
mass yields 4081 ± 117 signal events. With the to- 
tal reconstruction efficiency of 23.3%, the preliminary 
efficiency- corrected yield of D° — > n°n° based on our 
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i//(3770) data set is I752l±500(stat)±l559(syst) events 
( the estimation of the systematic uncertainty can be 
found in Ref. 1251 ). 




M k . [MeV/c 2 ) 

Figure 8: A fit to the D° — > n°n° candidate M/, c distribution based 
on the t/((3770) data sample. The black points are data, the black- 
smooth curve represents the overall fit (signal plus background), and 
the red-dashed curve corresponds to the fitted background shape. The 
blue-solid histogram represents the expected background shape and 
size based on our MC samples while the blue-dotted histogram is a fit 
to the data based on this expected MC-based background shape. 



We analyze D — > yy, and Fig.|9]shows AE distribu- 
tion based on the <//(3770) data set, where AE is the en- 
ergy difference between reconstructed energy of D me- 
son and beam energy. The signal candidates should be 
peak near zero, therefore, no significant signal events 
are observed. We perform a maximum-likelihood fit to 
the AE distribution. In the fit, the signal shape is fixed 
by the corresponding MC shape. The background shape 
consists of three parts; MC-based shape to represent the 
contamination from D° — > ttV whose size is also fixed 
based on our own observation; a 1st order polynomi- 
nal that covers the contamination from Bhabha events 
which appear smoothly over the entire AE spectrum; a 
1st order exponential polynominal, corresponding to the 
rest of the backgrounds. The fit yields -2.9 ±7.1 signal 
events. This translates into an upper limit of 1 1 events 
at 90% C.L. based on the Bayesian method. Including 
the estimated total systematic uncertainty, we arrive at 
S(D° -> yy)/S(D° -> ttV) < 5.8 x 1(T 3 at 90% C.L.. 
With the known value of S(D° -> 7T n°) 0, this corre- 
sponds to S(D° yy) < 4.7 x 10~ 6 . 

5. Conclusion 

Since the start of running in 2008, BESIII has taken 
about 2.9 ftr 1 of data at i/f(3770). With peak luminosity 




AE(MeV) 

Figure 9: A fit to the D° — > yy candidate AE distribution based on the 
i/<(3770) data sample. Black points are data, solid black curve is the 
overall fitted curve (signal plus backgrounds), the red-dashed curve is 
the fitted total background, and the green curve is the exponential and 
linear polynomials. 



reaching more than 6 x 10 32 cm _2 r' (60% of the de- 
signed luminosity), BESIII is poised to take more data 
at i/f(3770) and in the higher (D s ) energy region. BESIII 
collaboration presented the most precise measurement 
for D + — > yu + v decay. Using part of the data, BESIII 
has presented preliminary results of the D° — > K/nev 
decays. Results from the full dataset and other modes 
are coming in the near future. 

For the rare charm decay program, BESIII reported 
the search for D° — > yy decay with single tag technique, 
and backgrounds are expected to be lower than experi- 
ment at B factories. While we are waiting for BESIII 
to take more data at t/s = 3.773 GeV, there is an alter- 
nate analysis approach that is unique to our data sample. 
The produced i//(3770) in our sample decays into a pair 
of D°D°. Reconstructing one of the D° mesons with 
known exclusive modes while searching for D — > yy in 
the other D decay would yield an almost background- 
free environment, except for the irreducible contamina- 
tion from D° — > 7T°n° for which we have control. Such a 
study is also currently under way. We also explore other 
rare charm decays in the future with more dataset at the 
BESIII experiment. 
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